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Abstract 



The chiral phase transition at finite temperature is studied by using the 
Schwinger-Dyson equation in the dual Ginzburg-Landau theory, in which the 
dual Higgs mechanism plays an essential role on both the color confinement 
and the spontaneous chiral-symmetry breaking. At zero temperature, quark 
condensate is strongly correlated with the string tension, which is generated 
by QCD-monopole condensation, as (qq) 1 ^ ~ \fo. In order to solve the finite- 
temperature Schwinger-Dyson equation numerically, we provide a new ansatz 
for the quark self-energy in the imaginary-time formalism. The recovery of the 
chiral symmetry is found at high temperature; T c ~ lOOMeV with realistic 
parameters. We find also strong correlation between the critical temperature 
T c of the chiral symmetry restoration and the strength of the string tension. 
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In SU(iV c ) gauge theory, the appearance of magnetic monopoles was pointed out with the 
idea of the abelian gauge fixing based on the topological argument, 7r 2 (SU(A' c )/lJ(l) 7Vc " 1 ) = 
Z^ c_1 , by 't Hooft [[[]. He conjectured that the dual Meissner effect would be realized 
if QCD-monopoles were condensed [I]]. The basic idea for color confinement is that the 
color-electric flux between quarks is squeezed like a string or a tube, because the color- 
electric fields are excluded in the QCD-vacuum. This is similar to the Meissner effect in 
the superconductor, where the ordinary magnetic fields are excluded. The linear-confining 
potential is produced, since the squeezed flux has a uniform energy per unit length. The 
recent lattice QCD simulations support this idea that QCD-monopoles play a crucial role 
on color confinement through their condensation ||. 

The effective theory of non-perturbative QCD may be derived from QCD by using the 
abelian gauge fixing, which brings the non-abelian gauge theory into the abelian gauge the- 
ory with QCD-monopoles In order to describe the above mentioned dual Meissner effect, 
the dual Ginzburg-Landau (DGL) theory is constructed following Nambu's demonstration 
||, which uses the relativistic version of the Ginzburg-Landau theory on the Abrikosov vor- 
tex solution in the superconductor ||. QCD-monopole condensation causes strong and long 
range correlations between a quark and antiquark pair, which produce the linear-confining 
potential through the dual Higgs mechanism [|||| . We assume in addition the abelian dom- 
inance, in which the non-abelian part does not contribute to non-perturbative phenomena 
at low energy and is neglected J7|. We then studied chiral symmetry breaking in the DGL 
theory flJ^TOH . It was found that QCD-monopole condensation plays an essential role on 
the dynamical mass generation of the quark, which means spontaneous chiral-symmetry 
breaking, by using the Schwinger-Dyson (SD) equation for the dynamical quark |||8| |10 



We mention here that the lattice QCD study supports this idea of providing spontaneous 
chiral-symmetry breaking due to QCD-monopole condensation p"T| . |T2f . 

The high temperature and/or high density region would be realized in the laboratories 
as the intermediate states of the relativistic heavy-ion collisions, the study of which will be 
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made at RHIC in BNL and at LHC in CERN in the near future. The possible signals of 
a new phase, in which the liberation of colors, namely quarks and gluons, takes place, are 
expected to be seen in these experiments. Anticipating these experiments, the quark-gluon 
plasma (QGP) have been studied by many theorists. It is necessary for both confmement- 
deconfmement phase transition and chiral phase transition to understand the QCD phase 
transition. However, we do not have yet an effective theory to study the QGP physics 
based on the systematic treatment of these phase transitions, except for the lattice QCD 
simulations, which are not accessible for all the necessary physical quantities at the present. 

We study the QGP physics on the basis of the DGL theory, which describes both color 
confinement and chiral symmetry breaking at zero temperature [[|.|8] [L(J. In particular, we 
concentrate on the manifestation of chiral symmetry at finite temperature with confine- 
ment properties, which are based on the dual Meissner effect in the QCD vacuum. In this 
respect, we would like to mention the successful use of the DGL theory on the confinement- 



deconfinement phase transition at finite temperature |13|]. We shall formulate the finite- 
temperature SD equation using the imaginary-time formalism. 

The DGL theory is an infrared effective theory of non-perturbative QCD based on the 
dual Higgs mechanism in the abelian gauge . 

£ DGL = tr£ gauge (A M , + q(i? -eA- m q )q + tr[£> M , X ] f P>, x] - Atr( X f X - v 2 ) 2 , (1) 

where = d^ + igB^ is the dual covariant derivative. The dual gauge coupling g obeys the 
Dirac condition eg = An with e being the gauge coupling. The diagonal gluon and the 
dual gauge field B^ are defined on the Cartan subalgebra T 3 , T 8 on SU(3): = A^T 3 +AgT s , 
B^ = -B3T3 + BgT 8 . The QCD-monopole field x is defined on the SU(3) root vectors E a : 
X = XaE a - -S'gauge is the kinetic term of the gauge fields (A^, B^) in the Zwanziger 

form ||T4|| , 



^gaugc E--[n-(9A A)\ - ^[n • (d A B)} 2 — — [n • (d A A)\ v [n ■* (d A B)] v , (2) 

where the duality of the gauge theory is manifest. The DGL lagrangian has the [U(l) e ] 2 x 
[U(l) m ] 2 gauge symmetry. 



It is straightforward to calculate the static confining potential using the static quark- 
antiquark pair sources, which are totally color-singlet, located at a distance r = |b — a|. 

j M = Q^o{5 3 (x-b)-5 3 (x-a)} , (3) 

where Q = (Q 3 , Q 8 ) is the color charge of the static quark. The linear potential; or in the 
long distance is derived from the quark color-electric current-current correlation; —^i^D^ju 
with the static source (f|). Here, denotes the diagonal gluon propagator. We then get 



the string tension as the simple expression ||15|| , 



a - i , i ln ( — / A ) = Hi—* ln( — -rr^ (4) 

with Q 2 = • e 2 = ^ for N c = 3 and k = corresponding to the Ginzburg-Landau 



2tv c ^ 3 ^ ' c " " ^ V 3g 2 
parameter. The mass of B^, m# = V3gv, is proportional to the QCD-monopole condensate 

v For the type-II limit (n ^> V2), one finds a ~ 2irv 2 ln(2/€ 2 ), corresponding to the 

formula for the energy per unit length of the Abrikosov vortex in the type-II superconductor 
We may set parameters as e = 5.5, A = 25 and v = 0.126GeV, which reproduce 
the string tension {^/a ~ 0.44GeV) and the cylindrical radius of the hadron flux tube 
(R ~ m~ B l ~ 0.4fm) §. 

We get the ordinary SD equation by using the non-perturbative gluon propagator, ob- 
tained through QCD-monopole condensation in the DGL theory, as the full gluon propagator 
of QCD in order to include the non-perturbative effect in the infrared region. Hereafter, 
we concentrate on the case of the chiral limit for the simplicity of the argument. The SD 
equation for the dynamical quark propagator S q (p) in the rainbow approximation is written 
as 

s;\p) = ¥ + / -0^Q 2 i»s q (kh„D^( P - k) , (5) 

in the Euclidean metric. 

In the QCD-monopole condensed vacuum, the non-perturbative gluon propagator is de- 
rived by integrating out B^ fl,jg|-|T0|l , 



1 £\ l ial3£\vy6 n a n 'ykf3ks 



(6) 



k 2 k 2 + m 2 B (n-kf + a? 
where a e is the gauge fixing parameter on the residual abelian gauge symmetry [U(l) e ] 2 . It 
is noted that this modified gluon propagator is asymptotically reduced to the original gluon 
propagator of QCD in the ultraviolet region; k 2 ^> m 2 B . We have introduced the infrared 
cutoff a in this propagator (|5]) corresponding to the dynamical quark-antiquark pair creation 
and/or the color confinement in the size of hadrons [|10| . 

We take the angular average on the direction of the Dirac string n M in the SD equation, 
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(n • k) 2 + a? 



2vr 2 



1 



1 



(n ■ k)' 



a a + vF + 



a' 



(7) 



average 

Here, the dynamical quark is considered to move in various directions inside hadrons, and 
hence the constituent quark mass would be regarded as the quark self-energy in the angle- 



averaged case [ 10 1 . Taking a simple form for the quark propagator as S l (p) = if 1 — M(p 2 



the SD equation for the quark self-energy M(p 2 ) is obtained by taking the trace of Eq.( 
in the Landau gauge (a e = 0), 



M(p 2 ) 



d A k 



x 



Q 2 



M(k 2 



(2tt) 4 ^ k 2 + M 2 (k 2 
2 



k 2 + m 2 B k 2 



1 4 
+ — + 

a a + \/k 2 + a 2 



' m 2 B — a 2 a 2N 



k 2 + 



B 



k 2 



(8) 



with kfj, = — k^. It is noted that the r.h.s. of Eq.(H) is always non-negative. 

We solve Eq.(^|) numerically using the Higashijima-Miransky approximation [16], Q 2 = 
-i-nCp ■ a s eff (max{p 2 , k 2 }) with Cp = f° r SU(A^ C ) in order that the SD equation is 

reduced to the usual one of the QCD-like theory |fT6lJT7| in the ultraviolet region. As a 



consequence, the quark self-energy M(p 2 ) has the asymptotic form, which is obtained by a 
study of the operator product expansion and the renomalization group [IB], in the ultraviolet 
limit. The running coupling with the hybrid type behavior is obtained as 



a s ett (p 2 ) 



12tt 



A 



48tt 2 (A c + 1) 



QCDeXp[ (llA c -2A / )-e 2 



(9) 



(llN c -2N f )ln[(p 2 +p 2 c )/A 2 lCD ] ' ^v c -^ v/ 
where p c approximately divides the momentum scale into the infrared region and the ultra- 
violet region with the QCD scale parameter Aqcd fixed at 200MeV. Here, p c is fixed by 



a s eS (p 2 = 0) = a s j (N c + 1); a s = e 2 /47r. The magnitude of the quark condensate plays the 
role of an order parameter of spontaneous chiral-symmetry breaking. We can calculate the 
quark condensate with the thus obtained quark self-energy M(p 2 ) from the SD equation as 

/ ^ N c f A2 „ 2 k 2 M{k 2 ) 

where the ultraviolet cutoff A is introduced to regularize the ultraviolet divergence of the 
integral. If the ultraviolet cutoff is taken in the deeply asymptotic region, A ^> tub, one 
can get the quark condensate in the renormalizat ion-group invariant form for N c = Nf = 3, 
{qq)rgi = (QQ) A ' {^(AVAqcd)} The ultraviolet cutoff is taken to be large enough 
as A/Aqcd = 10 3 hereafter. One finds that QCD-monopole condensation provides a large 
contribution to spontaneous chiral-symmetry breaking |]6|j8|-[l0| , because the mass of the 
dual gauge field, rriB, is proportional to the QCD-monopole condensate |5]|| and the quark 
condensate increases with m# [|G|,|3|-|T0|| . In other words, the strength of spontaneous chiral- 
symmetry breaking seems originated from QCD-monopole condensation by the sufficient 
strength of the string tension as shown in Fig. 1. The parameter set is taken as e = 5.5, 
ins — 0.5GeV and a = 85MeV, where these parameters except for a are set by fitting the 
confinement phenomena at zero temperature. In this case, the quark condensate and the 
pion decay constant at zero temperature are well-reproduced as (qq) RGl — — (247MeV) 3 and 
f n ~ 88MeV, respectively |0 . 



Now we study the chiral symmetry restoration at finite temperature using the imaginary- 
time formalism for the SD equation [Q. The finite-temperature SD equation is obtained 
by making the following replacement in Eq. (H) [[0|] : 



P4 — > oJn = (2n + 1)ttT , (Ha) 
f d 4 k « r dk 

/ 7^4 ^ T L / 7^3 > ( llb ) 
(27T) n=-oo J (27T) 

M(p 2 ) ^ M T (u n ,p) . (11c) 

The resulting equation is very hard to solve even numerically since the quark self-energy 

depends not only on the three dimensional momentum p, but also on the Matsubara fre- 



quencies uo n ||20|| . We propose the covariant-like ansatz [181, in which a replacement is made 
for the quark self-energy at T 7^ instead of ( |llc| ) as 

M(p 2 ) - M T (p 2 + ^) ^ M T (p 2 ) (12) 



with p 2 = p 2 + u; 2 and ui n = (2n + 1)ttT fll8 |. It is noted that this ansatz guarantees that the 
finite-temperature SD equation in the limit T — > is exactly reduced to the SD equation 
(§) at T = 0. This fact is also confirmed by numerical calculations as shown in Fig. (2). The 
final form of the SD equation at T 7^ is derived as 



M T (p 2 



— y 



IT*' ' ' I 2 
" m=— 00 ^ 



00 A pi 
2 



M T (k 2 



k 2 + M 2 (k 2 ) 



2 1 4 

fcH™ + ?77,r & 2 a 



in 



B 



+ \/ + a 2 \ knm + m l ^nr 



(13) 



where k 2 m = k + p 2 — 2z\j (p 2 — uj 2 )(k 2 — uj 2 ^) — 2uj m u n . It is noted that the large number 
of the Matsubara frequencies have negligible contribution at high temperature in Eq.(pT3[) 



since k 2 



~ (UJ m — LO n 



at T ^> |p|, |k|. In other words, it is important for the small number 



of the Matsubara frequencies to generate the non-trivial solution of the quark self-energy 
M T {p 2 ) in r.h.s. of Eq. ([l3|) . We then solve the SD equation by setting u n = in the r.h.s. 
of Eq . (|13"D . In addition, it is necessary to truncate the infinite sum of m in order to solve it 
numerically. We check the numerical convergence by this truncation for the m sum, which 
is quite fast at high temperatures. In the effective running coupling eq.(|9|), there is only a 
slight modification as Q 2 = AtiCf • af 1 (ma.x{p 2 , k 2 }) due to the covariant-like ansatz. We 
show in Fig. 2, the quark self-energy M T (p 2 ) as a function of p 2 at finite temperature. No 
nontrivial solution is found in the high temperature region, T ~ HOMeV. In other words, 
the chiral symmetry is restored at high temperature. 

The quark condensate is easily calculated using the quark self-energy M T (k 2 ) as 



^raax p 

—AN C ■ T Y, J 
2N C 



a dk 



M T (u; n ,k) 



(27ryk 2 + u 2 n + M 2 (uj n ,k) 



7T Z 



dk 2 \ k 2 -LU 2 - 



M T {k 2 



n=Q ""n 



'k 2 + M 2 (k 2 ) 



(14) 



with n max = — |]. The quark condensate, {qq) T , shown in Fig. 3, decreases gradually 
with temperature in the low temperature region (~ 0.5T C ) and vanishes suddenly near the 
critical temperature T a . We mention that the behavior of the quark condensate in the low 
temperature region is similar to that of the calculations on the chiral perturbation theory, 
namely (qq) T / (qq) T=0 = 1 — T 2 /8f% + 0(T 4 ) f2lf . In addition, these results are slightly 
different from the lattice QCD simulation p2[, which shows that the temperature dependence 



of the chiral condensate is small until near the critical temperature. In these calculations, 
we do not take any temperature dependence in the parameters used for simplicity, although 
we should take it into account as more information on high temperature QCD becomes 
available. 

Finally, we examine the correlation between the critical temperature T c of the chiral 
symmetry restoration and the confinement quantity such as the string tension a. As shown 
in Fig. 4, there exists a strong correlation between them. From the physical point of view, 
higher T c would be necessary to restore the chiral symmetry in the system where the linear- 
confining potential is much stronger, because there exists the close relation P,|l0|-|T2|1 between 
color confinement and dynamical chiral-symmetry breaking through QCD-monopole conden- 
sation as shown in Fig. 1. The critical temperature is estimated as lOOMeV ~ T c ~ HOMeV, 
when one takes the standard value of the string tension (yfa ~ 0.44GeV). We remark here 
that the critical temperature would increase if we take smaller cutoff parameter a at finite 
temperature, which is expected as due to the increase of the size of hadrons because of 
smaller confining force |13|]. Hence, it would be very interesting to couple the calculations 
of QCD-monopole condensation and chiral symmetry breaking at finite temperature. 

In summary, we have studied spontaneous chiral-symmetry breaking and its restoration 
at finite temperature using the SD equation in the DGL theory, which is known to provide 
both quark confinement and dynamical chiral-symmetry breaking. We have solved the finite- 
temperature SD equation numerically with the covariant-like ansatz. The chiral symmetry 
is restored at T c ~ lOOMeV. We have found the strong correlation between the critical 
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temperature T c and the string tension a. 

We acknowledge fruitful discussions with O. Miyamura, S. Umisedo and H. Ichie on the 
relation between spontaneous chiral-symmetry breaking and monopole condensation. 
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FIGURE CAPTIONS 



FIG.l. The correlation between the quark condensate and the mass of the dual 
gauge field and also the string tension. The other parameters are fixed 
as e = 5.5 and a = 85MeV. The unit is taken by A QC d - 200MeV. 

FIG. 2. The quark self-energy M T (p 2 ) as a function of p 2 at T = 0, 60 and 
lOOMeV with e = 5.5, m B = 0.5GeV and a = 85MeV. 

FIG. 3. The ratio of the quark condensate at finite temperature to that of zero 
temperature as a function of temperature normalized to the critical tem- 
perature T c . The same parameters are used as in Fig. 2. T c is found at 
T c ~ HOMeV in this case. 

FIG. 4. The critical temperature T c as a function of the square root of the string 
tension a. T c rises almost linearly with yfd. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9606313vl 



This figure "figl-2.png" is available in "png" format from: 
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